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ABSTRACT 

We present medium resolution spectroscopy and multi-epoch VRI photometry for 21 new nearby 
(< 50 pc) white dwarf systems brighter than V ^ 17. Of the new systems, ten are DA (including a 
wide double degenerate system with two DA components), eight are DC, two are DZ, and one is DB. 
In addition, we include multi-epoch VRI photometry for eleven known white dwarf systems that do 
not have trigonometric parallax determinations. Using model atmospheres relevant for various types 
of white dwarfs (depending on spectral signatures), we perform spectral energy distribution modeling 
by combining the optical photometry with the near-infrared JHKs from the Two Micron All-Sky 
Survey to derive physical parameters (i.e., Toff and distance estimates). We find that twelve new 
and six known white dwarf systems are estimated to be within the NStars and Catalog of Nearby 
Stars horizons of 25 pc. Coupled with identical analyses of the 56 white dwarf systems presented in 
Paper XIX of this series, a total of 20 new white dwarf systems and 18 known white dwarf systems 
are estimated to be within 25 pc. These 38 systems of the 88 total studied represent a potential 
34% increase in the 25 pc white dwarf population (currently known to consist of 110 systems with 
trigonometric parallaxes of varying qualities). We continue an ongoing effort via CTIOPI to measure 
trigonometric parallaxes for the systems estimated to be within 25 pc to confirm proximity and further 
fill the incompleteness gap in the local white dwarf population. Another 38 systems (both new and 
known) are estimated to be between 25 and 50 pc and are viable candidates for ground-based parallax 
efforts wishing to broaden the horizon of interest. 

Subject headings: solar neighborhood — stars: distances — stars: evolution — stars: statistics — 
white dwarfs 



1. INTRODUCTION 

Given that all stars less than '^8 Mq will eventually 
become white dwarfs (WDs), the study of this class of 
objects is vital to understanding our Galaxy. In particu- 
lar, WD research addresses questions concerning stellar 
structure and evolution. Galactic components (i.e., thin 
disk, thick disk, halo), and even dark matter. The old- 
est (i.e., coolest and least luminous) WDs help to con- 
strain the age of the Galactic components (particularly 
the thick disk) but shine feebly and are only visible and 
easily characterized if they are in our solar neighbor- 
hood. A complete volume-limited sample provides ac- 
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curate statistics that can be reliably applied to the rest 
of the Galaxy to infer the WD mass fraction, contribu- 
tion to the halo population, and WD number density 
within the disk. Interesting candidates can be targeted 
and more thoroughly scrutinized to identify unusual and 
astrophysically compelling systems. 

In a continuing effort to further complete the nearby 
WD sample, we present spectra, optical and near- 
infrared photometry, as well as modeled physical param- 
eters for 21 new WD systems in the southern hemisphere 
brighter than V ~ 17. Of these, 12 are estimated to be 
within 25 pc, the horizon of the Catalog of Nearby Stars 
(CNS, Ghese & Jahreifi 1991) and the NStars Database 
(Henry et al. 2003). In addition, six previously known 
WDs without trigonometric parallaxes are estimated to 
be within 25 pc. 

2. CANDIDATE SELECTION 

We conducted a southern hemisphere proper motion 
search using digitized scans of the SuperCOSMOS Sky 
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Fig. 1. — Reduced proper motion diagram used to select WD 
candidates for spectroscopic follow-up. Plotted are the 306 new 
high proper motion objects from Subasavage et al. (2005a,b) for 
which both i?59F and J magnitudes were available. The line is a 
somewhat arbitrary boundary between the WDs (below) and the 
subdwarfs (just above). Main sequence dwarfs fall above and to 
the right of the subdwarfs, although there is significant overlap. 
Asterisks indicate the 23 new WDs reported here. The six filled 
triangles in the WD region are confirmed WDs from the SCR sur- 
vey and were published in Subasavage et al. (2007). The encircled 
point labeled "sd" is a confirmed subdwarf contaminant in the WD 
region. 

Survey (SSS), adopting a faint magnitude limit of R^gp = 
16.5, ealled the SuperCOSMOS-RECONS (SCR) survey. 
The first wave of the survey adopted a lower proper mo- 
tion limit of 0.40" yr~^ with new discoveries published 
in Hambly ct al. (2004), Henry ct al. (2004), and Sub- 
asavage ct al. (2005a, b). A second effort conducted by 
Finch et al. (2007) surveyed objects with proper motions 
from 0.18-0.40" yr~^ and covered declinations —47° to 
—90°. The WD candidate selection effort also adopted 
the lower proper motion limit of 0.18" yr^^, yet reached 
farther north to declination = 00°. Thus, this effort cov- 
ered 92% of the southern sky, avoiding a few regions 
near the Galactic plane and the Magellanic Clouds (see 
Figure 1 of Subasavage et al. 2005b). A full discus- 
sion of the search methodology can be found in Ham- 
bly et al. (2004). Briefly, each of the four plate scans 
used (Bj. i?ESOi -R59FJ ^ivn) were positionally mapped to 
a common coordinate system. Any object that appeared 
on all four plates and had a proper motion less than the 
lower limit was discarded. Objects that remained were 
then searched out to a radius defined by a proper motion 
of 10.00" yr~^ and over 360° for any other unpaired ob- 
jects. After a series of automated sifts for false pairings, 
the remaining objects were cross-referenced with previ- 
ous proper motion catalogs to identify new and known 
objects. 

Near-infrared JHKs magnitudes were extracted from 
2MASS for all real objects (new and known) and a re- 
duced proper motion (RPM) diagram using the i^sgp — J 
color was generated (see Figure 1). RPM (in this case, 
designated by Hr^^^ because the iJsgp magnitude was 
used) is a quantity similar to absolute magnitude and is 
defined by 

Hr,„^ = i?59F + 5 + 5 log(/x) (1) 

where /i is proper motion. It serves to relate two ob- 
served quantities, apparent magnitude and proper mo- 
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Fig. 2. — Plot of infrared J—Ks color (transformed from the CIT 
system to the 2MASS system using the transformation equations 
of Carpenter 2001) vs. My for WDs within 25 pc (Bergeron et 
al. 2001) and late-type subdwarfs within 60 pc (W.-C. Jao 2008, 
private communication) . Note how this particular color is less than 
0.5 {vertical dashed line) over all My and is degenerate for the 
WDs. Numbered points are discussed in the text. 

tion, with two intrinsic quantities, luminosity (absolute 
magnitude) and tangential velocity. The advantages of 
using the J?59f — J color is that all SCR detections have 
an R5QF magnitude (defined by the survey). Given that 
WDs are relatively blue, the likelihood of registering a 
near-infrared magnitude in the 2MASS database is great- 
est at J because its limiting magnitude is ~1.0 mags 
fainter than Kg. Also, this color incorporates both op- 
tical plate and near-infrared magnitudes that minimize 
the intrinsic uncertainties with plate magnitudes (e.g., 
non-linearity) . Both new SCR discoveries and known re- 
coveries fell within the WD region of the RPM diagram. 
Most of the known objects had already been classified as 
WDs by previous researchers. Yet, a significant number 
of known high proper motion (HPM) objects within this 
region were not classified as WDs. It is likely that they 
escaped identification because of the poor plate magni- 
tudes previously available. For instance, all plate magni- 
tudes fainter than rripg = 10.0 in the classic Luyten Half- 
Second (LHS) Catalogue (Luyten 1979a) and the New 
Luyten Two-Tenths (NLTT) Catalogue (Luyten 1979b) 
are by-eye estimates determined by Luyten (1979a). Rig- 
orous calibrations performed by the SSS team yield high 
quality plate magnitudes that reveal a number of new, 
relatively bright, WDs. 

As an additional constraint, especially for objects near 
the arbitrary boundary that delineates the subdwarfs 
and the WDs, a J — Kg color sift was implemented. As 
is evident in Figure 2, this color is degenerate for WDs 
(i.e., there is no unique absolute magnitude for a given 
color) and is always less than 0.5 mag (for single WDs). 
In contrast, subdwarfs near the boundary typically have 
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J — Ks ^ 0.6 or larger and thus there exists a signifi- 
cant color gap to distinguish between the two luminosity 
classes. Indeed, four new WDs were spectroscopically 
identified that are found in the subdwarf region of the 
RPM diagram. Two exceptional points in the WD region 
are: (1) WD 0548-001, which is a DQ WD (discussed 
in Dufour et al. 2005) that has a magnetic field of ~10 
MG (Kawka et al. 2007); (2) WD 2251-070, which is a 
DZ WD (discussed in Dufour et al. 2007) that is heavily 
line blanketed (Wesemael et al. 1993, Figure 11). 

Before the candidates were targeted for spectroscopic 
follow-up observations, the basic linear plate relation of 
Oppenheimer et al. (2001) was used to estimate dis- 
tances to candidates, assuming they were WDs. Only 
those candidates whose distance estimates were within 
~30 pc were selected for spectroscopic observations (this 
step was omitted from the procedures resulting in the 
new WD discoveries presented in Subasavage et al. 2007, 
hereafter referred to as Paper XIX). This constraint fa- 
vored cooler, nearby stars so that even though the sample 
size of new systems presented here is smaller than that 
of Paper XIX (21 vs. 33), more systems are estimated 
to be within 25 pc (twelve vs. eight). A total of 21 new 
WD systems (containing 22 WDs) were spectroscopically 
confirmed and are hereafter referred to as the "new sam- 
ple" . Optical photometry observations were obtained to 
improve distance estimates (discussed in § 4.1). 

Eleven known WDs without trigonometric parallaxes 
(hereafter referred to as the "known sample") were also 
targeted for optical photometry observations to improve 
their distance estimates. The known sample includes 
nine objects that previous authors predicted to be nearby 
(i.e., Aannestad et al. 1993: Holberg et al. 2002, 2008: 
Kawka et al. 2004, 2007; Pauh et al. 2006; Venues & 
Kawka 2003) as well as two objects found via RPM that 
were known WDs but whose distances had not yet been 
estimated. 

3. DATA AND OBSERVATIONS 

3.1. Astrometry and Nom.enclaf.ure 

In keeping with the traditional naming convention for 
WDs, which uses the object's epoch 1950 equinox 1950 
coordinates, 2MASS coordinates (when available, other- 
wise SSS coordinates were used) for the new systems pre- 
sented here were extracted and adjusted for proper mo- 
tion from the epoch of observation to epoch 2000 equinox 
2000. The coordinates were then precessed to equinox 
1950 using the IRAF task precess and again adjusted for 
proper motion (opposite in direction) to yield epoch 1950 
equinox 1950 coordinates. 

Proper motions for both the new and known samples 
were taken primarily from the SCR proper motion sur- 
vey, but in a few cases from elsewhere in the literature. 
The Appendix contains the proper motions used for co- 
ordinate adjustment, as well as J2000.0 coordinates and 
alternate names. 

3.2. Spectroscopy 

Spectroscopic observations were taken during two runs 
in May and December 2006 at the Cerro Tololo Inter- 
American Observatory (CTIO) 1.5 m telescope as part 
of the Small and Moderate Aperture Research Telescope 
System (SMARTS) Consortium. The Ritchey-Chretien 
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Fig. 3.— Spectral plots of the hot (Teff > 10000 K) DA WDs 
from the new sample, plotted in descending Toff as derived from the 
SED fits to the photometry. Spectra have been normalized at 5200 
A and offset by integer amounts. The inclusion of the spectrum of 
known object WD 2133—135 in this plot is discussed in § 4.2. 

spectrograph and Loral 1200 x 800 CCD detecter were 
used with grating 09 (in first order), providing 8.6 A res- 
olution and wavelength coverage from 3500 to 6900 A. A 
slit width of 6" was used and prevented the preferential 
light loss, at either the blue or the red end, encountered 
in the data presented in Paper XIX. Observations con- 
sisted of two exposures (typically 20-30 minutes each) 
to permit cosmic ray rejection, followed by a compari- 
son HcAr lamp exposure to calibrate wavelength for each 
object. Bias subtraction, dome/sky flat-fielding, and ex- 
traction of spectra were performed using standard IRAF 
packages. 

Spectroscopic flux standards (two or three) were ob- 
served each night to calibrate the response of the CCD 
across the dispersion axis. However, because the selected 
flux standards are bright, neutral density filters of either 
2.5 or 5.0 magnitudes extinction were used. Thus, the re- 
sulting science spectra are relatively flux calibrated (i.e., 
the slopes are correct) but not on an absolute flux scale 
(i.e., erg cm~^ s~^ A~^). All spectra are normalized at 
5200 A for the sake of plotting. 

Spectra for the new DA WDs with Tcff > 10,000 K 
are plotted in Figure 3, while spectra for new DA WDs 
with Teff < 10,000 K are plotted in Figure 4. Spectra 
for the new featureless DC WDs are plotted in Figure 5. 
Spectra, as well as model fits, for two new calcium-rich 
DZs, the hottest DA, and a helium-rich DB are shown in 
Figures 6 and 7 and are described in § 4.1. 

3.3. Photometry 

Optical Johnson-Kron-Cousins VRI^ photometry for 
the new and known WDs was obtained at the CTIO 
0.9 m telescope during several observing runs from Jan- 
uary 2006 through March 2008 as part of the SMARTS 
Consortium. The 2048 x 2046 Tektronix CCD cam- 
era was used with the Tex 2 VRI filter set. Standard 
stars from Graham (1982) and Landolt (1992, 2007) were 
taken nightly through a range of airmasses to calibrate 
fluxes to the Johnson-Kron-Cousins system and to cal- 
culate extinction corrections. Bias subtraction and flat- 
fielding (using calibration frames taken nightly) were per- 
formed using standard IRAF packages. When possible, 

^ The central wavelengths for Vj, Rkc, and /kg are 5475, 6425, 
and 8075 A, respectively. 
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Fig. 4.— Spectral plots of cool {T^s < 10000 K) DA WDs 
from the new sample, plotted in descending T^g as derived from 
the SED fits to the photometry. Spectra have been normalized at 
5200 A and offset by integer amounts. The asymmetry in the Ha 
feature in the spectrum of WD 2009—471 is the result of a cosmic 
ray landing just redward of Ha that could not be reliably removed. 

an aperture of 14" in diameter (consistent with Landolt 
1992) was used to determine stellar flux. If cosmic rays 
fell within this aperture, they were removed before flux 
extraction. In cases of crowded fields or nearby com- 
panions, aperture corrections were applied and ranged 
from 4" to 12" in diameter using the largest aperture 
possible without including contamination from neighbor- 
ing sources. Uncertainties in the optical photometry are 
±0.03 mag in each band and incorporate both internal 
(night-to-night variations) and external (fits to the stan- 
dard stars) uncertainties^. The final optical magnitudes 
are listed in Table 1, as well as the number of nights each 
object was observed. 

4. ANALYSIS 

4.1. Modeling of Physical Parameters 

A detailed description of the model atmospheres used 
to model the WDs can be found in Paper XIX and refer- 
ences within. Briefly, the optical VRI and 2MASS near- 
infrared JHKs magnitudes are converted into observed 
fluxes and compared to the resulting SEDs predicted by 
the model atmospheres calculations. The observed flux, 
/™, is related to the model flux by the equation 



/r = 47r {R/D f HI 



(2) 



where R/D is the ratio of the radius of the star to its dis- 
tance from Earth, and if™ is the Eddington flux (depen- 
dent on Toff, log 5, and atmospheric composition), prop- 
erly averaged over the corresponding filter bandpass. We 

^ A complete discussion of the error analysis can be found in 
Henry et al. (2004). 



taO 



T3 
(D 
N 



O 

S5 



1116-470 
0431 —279 

' ^11 Q_TQO 



0751-252 



Wavelength (A) 
Fig. 5. — Spectral plots of featureless DC WDs from the new 
sample, plotted in descending T^g as derived from the SED fits 
to the photometry. Spectra have been normalized at 5200 A and 
offset by integer amounts. 

consider only Tes and the solid angle [tt{R/D)'^] to be free 
parameters, and the uncertainties of both parameters are 
obtained directly from the covariance matrix of the fit. 
In this study, we simply assume a value of log 5 = 8.0 for 
each star. As discussed in Bergeron et al. (1997, 2001), 
the main atmospheric constituent — hydrogen or helium 
— is determined by comparing the fits obtained with 
both compositions, or by the presence of Ha in the op- 
tical spectra. The derived values for T^g for each object 
are listed in Table 1. For the two new DZ stars, spectro- 
scopic modeling of metal lines as well as the photometry 
was used to constrain Tes (see below). Also listed are the 
spectral types for each object determined based on their 
spectral features. The DAs have been assigned a half- 
integer temperature index as defined by McCook & Sion 
(1999), where the temperature index equals 50,400/Teff. 

Our grids of model atmospheres and synthetic spectra 
for DA and DB stars are described respectively in Licbert 
et al. (2005) and Beauchamp et al. (1996). The atmo- 
spheric parameters - Toff, log g (and hydrogen abundance 
for the DB stars) - are determined from the optical spec- 
tra using the so-called spectroscopic technique (see, e.g., 
Bergeron et al. 1992), which relies on a detailed com- 
parison between synthetic and observed normalized line 
profiles. The model spectra (convolved with a Gaussian 
instrumental profile) and the optical spectrum are first 
normalized to a continuum set to unity. The calculation 
of is then carried out in terms of these normalized 
line profiles only, and the best fitting solution is obtained 
using the nonlinear least-squares method of Leveiibc;rg- 
Marquardt (Press et al. 1992), which is based on a steep- 
est descent method. Spectroscopically derived values of 
Toff and log g are given in the footnotes of Table 1 for 
the DA and DB stars. 

For the DZ stars' spectra in Figure 6, we rely on the 
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TABLE 1 

Optical and Infrared Photometry, and Derived Parameters for New and Known White Dwarfs. 



No. 
of Obs. 



J eft 

(K) Comp. 



(PC) 



Spec. 
Type 



New Spectroscopically Confirmed White Dwarfs 



0011-721 . . . 
0102-579 . . . 

0123-460 . . . 
0134-177. . . 
0149-723 . . . 
0311-649 . . . 
0431-360 . . . 
0431-279 . . . 
0620-402 . . . 
0651-398A . 
0651-398B . 
0655-390 . . . 

0708- 670 . . . 

0709- 252 . . . 
0751-252 . . . 
0816-310 . . . 
0856-007 . . . 
1116 — 470 . . . 
1817-598 . . . 
1916-362 . . . 
2009-471 . . . 
2118-388 . . . 
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^Too cool for a reliable spectroscopic fit (i.e., minimal Balmer line absorption). 

^Spectral fit yielded a Tgff = 8228 ± 138 K and log g = 8.19 ± 0.12. Common proper motion companion to NLTT 3566 (see § 4.2). 

•^Spectral fit yielded a T^.^.f = 11613 ± 192 K and log g = 8.14 ± 0.06. 

''spectral tit yielded a T^g = 12632 ± 300 K and log g = 7.63 ± 0.06 (see § 4.2). 

"^Spectral fit yielded a T^ff = 7214 ± 135 K and log g ^ 7.68 ± 0.19. Common proper motion companion to WD 0651-398B (see § 4.2). 
^Spectral fit yielded a T^ff = 11208 ± 180 K and log g - 8.12 ± 0.06. 
^Common proper motion companion to LTT 2976 (see § 4.2). 

'^Spectral fit yielded a Tgff = 27830 ± 1050 K and log g = 7.97 ± 0.07 (see § 4.2 and Figure 7). 

Vennes &^ Kawka (2003) estimate a distance of 15 pc. 
^Holberg et al. (2002) estimate a distance of 10.79 pc. Holberg et al. (2008) estimate a distance of 12.05 pc. 

Holberg et al. (2002) estimate a distance of 18.22 pc. 

'spectral fit yielded a Tg££ = 10357 ± 158 K and log g = 7.99 ± 0.07. Pauli et al. (2006) estimate a distance of 23.3 pc. 
Kawka et al. (2007) estimate a distance of 14 pc. Holberg et al. (2008) estimate a distance of 14.24 pc. 
^Aannestad et al. (1993) estimate a distance of 28 pc. 
*^Kawka et al. (2004) estimate a distance of 33 pc. 

^Holberg et al. (2002) estimate a distance of 15.6 pc. Holberg et al. (2008) estimate a distance of 17.45 pc. 
'^Holberg et al. (2002) estimate a distance of 12.41 pc. 



procedure outhned in Dufour et al. (2007). We obtain 
a first estimate of the atmospheric parameters by fitting 
the photometric SED with an assumed value of the metal 
abundances (assuming solar abundance ratios). We then 
fit the optical spectrum to measure the metal abun- 
dances, and use these values to improve our atmospheric 
parameters from the photometric SED. This procedure 
is iterated until a self-consistent photometric and spec- 
troscopic solution is achieved. In the cases of the two DZ 
stars presented here, a slight abundance of hydrogen was 
included [log (H/He) = — 3] in the modeled fits. As shown 
in Dufour et al. (2007), the Ca line widths and depths 
are affected by the presence of hydrogen even if the abun- 
dance is not significant enough to produce spectral sig- 
natures. Our spectra were best fit with the inclusion of 
hydrogen. The values of T^a listed in Table 1 for the new 
DZ stars were derived using this additional spectroscopic 



modeling constraint [i.e., "He(-f Ca)" in the composition 
column]; whereas, the known DZ (WD 2216—657) was 
modeled using only photometry assuming a pure He at- 
mosphere (i.e., "He" in the composition column). 

Once the effective temperature and the atmospheric 
composition were determined, we calculated the abso- 
lute visual magnitude using a procedure identical to that 
outlined in Paper XIX (i.e., combined the new photomet- 
ric calibration of Holberg et al. 2006, with evolutionary 
models^). We then compared the absolute magnitude to 
the apparent V magnitude observed to derive a distance 
estimate for each star (reported in Table 1). Errors on 
the distance estimates incorporate the errors of the pho- 
tometry values as well as an error of 0.25 dex in log g, 
which is the measured dispersion of the observed distri- 

* See http: / /www.astro.umontreal.ca/~bergeron/CoolingModels/. 
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Fig. 6. — Spectral plots of the two DZ WDs from the new 
sample. The inset plots display the spectra {thin lines) in the re- 
gions to which the models {thick lines) were fit. Spectroscopic best 
fit physical parameters are listed below. Prominent lines from Ca 
{vertical lines), Mg {open circles), and Fe {filled circles) arc labeled. 
Both modeled fits incorporated a slight hydrogen abundance of log 
(H/He) = -3 (see § 4.1). 

bution using spectroscopic determinations (see Figure 9 
of Bergeron et al. 1992). 

4.2. Comments on Individual Systems 

WD 0102-579: A new HPM object discovered dur- 
ing the SCR proper motion survey (/i = 0.239" yr~-^ at 
position angle 91.1°, Finch et al. 2007). It is likely 
a common proper motion companion to a previously 
known HPM object, the red dwarf NLTT 3566 {fi = 
0.257" yr-i at position angle 87.5°), separated by 110.6" 
at position angle 206.6°. The Hipparcos parallax for 
NLTT 3566 is 19.51 ± 3.13 mas (distance = 51.3 ± 9.8 
pc, van Leeuwen 2007) - entirely consistent with our dis- 
tance estimate for the WD of 44.4 ± 7.5 pc. 

WD 0311-649: The hottest DA WD in the new sam- 
ple (Teff = 11,945 ± 557 K). The Balmer lines of our high 
S/N spectrum are best fit to a spectral fitting model with 
Teff = 12,632 ± 300 K and log g = 7.63 ± 0.06 (see Figure 
7), significantly less than the assumed value of log g = 8.0 
in the photometric analysis. This would imply that the 
object is low mass (M = 0.42 ± 0.03 Mq) and more lu- 
minous thereby making it more distant (~28 pc). Given 
the age of our Galaxy, the lowest mass WD that could 
have formed is -0.47 Afg (Iben & Renzini 1984). If the 
mass is correct, it is extremely unlikely that this object 
formed through single-star evolution and thus is likely a 
multiple system. Trigonometric parallax measurements 
are underway to confirm the luminosity and hence the 
mass via our Cerro Tololo Inter- American Observatory 
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Fig. 7. — Top: Spectral plot of the DA (hydrogen-rich) WD 
0311—649 with the model fit {thick line) overplotted. Spectro- 
scopic best fit physical parameters are listed below. Bottom: Spec- 
tral plot of the DB (helium-rich) WD 1916-362. The inset plot 
displays the spectrum {thin line) in the region to which the model 
{thick line) was fit. Spectroscopic best fit physical parameters are 
listed below. Note that the sharp absorption feature at ~4100 A is 
an artifact produced by a cosmic ray that could not be reliably 
removed because of its proximity to a true spectral feature. 

Parallax Investigation (CTIOPI, Jao et al. 2005; Costa 
et al. 2005, 2006; Henry et al. 2006). 

WD 0620-402: A new DZ WD whose modeled phys- 
ical parameters place this object in a lower temperature 
realm where high atmospheric pressure effects exist that 
are not included in the Ca-rich (DZ) models described in 
§ 4.1. Thus, the model fit is likely inaccurate for this ob- 
ject, so that the physical parameters derived are not well 
constrained (especially given the low S/N of the spec- 
trum, see Figure 6). 

WD 0651-398AB: A widely separated (87.2" at po- 
sition angle 317.4°) HPM double degenerate system dis- 
covered during the SCR proper motion survey. The A 
component {n = 0.229" yr~^ at position angle 344.7°) 
has Toff = 7222 ± 219 K and the B component = 
0.227" yr^i at position angle 344.4°) has T^s = 6450 ± 
220 K. A third object, WT 204, is separated by 59.3' 
at position angle 67.2° and has a similar proper motion 
of 0.213" yr-i at position angle 341.9°. We have ob- 
tained spectra and optical photometry, which indicate a 
spectral type of M3.0V and Vj = 13.17, Rkc = 12.03, 
Ikc = 10.60. Using the main-sequence distance relations 
of Henry et al. (2004), we obtain a distance estimate of 
19.8 ±3.1 pc, which is entirely consistent with the dis- 
tance estimates to the A and B components, 25.1 ± 4.3 
pc and 26.9 ± 4.5 pc, respectively. While highly unlikely 
that this is a bound system (the third component's sep- 
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aration is ^ 90,000 AU), it is possible that they are part 
of a moving group. 

WD 0751-252: A new HPM object discovered dur- 
ing the SCR proper motion survey (/z = 0.426" yr~^ at 
position angle 300.2°, Subasavage et al. 2005b). Its 
proper motion is similar to a previously known HPM 
object, LTT 2976 {fi = 0.361" yr^^ at position angle 
303.8°). The distance estimate for the WD (17.8 ± 2.9 
pc) is consistent with the Hipparcos parallax of 51.53 ± 
1.46 mas (distance = 19.41 ± 0.57 pc, van Leeuwen 2007) 
for LTT 2976. Thus, these two objects likely form a sys- 
tem with a separation of 6.6' (~ 8000 AU) at position 
angle 208.9°. 

WD 0816-310: A new DZ WD whose spectrum has 
been reliably reproduced using the methods appropriate 
for DZs described in § 4.1 (see Figure 6). The inset plot 
extends down to 3700 A (rather than 3800 A for all other 
spectral plots) to illustrate the validity of assuming solar 
abundance ratios (at least for Ca, Mg, and Fe) given the 
quality of the fit. 

WD 1916-362: The only DB WD in the new sam- 
ple. The photometric temperature (Toff = 24,105 ± 8797 
K) agrees reasonably well with the spectroscopic temper- 
ature (Teff = 27,830 ± 1050 K) assuming a pure helium 
composition given that no hydrogen features are visible 
in the spectrum. At these temperatures, trace amounts 
of hydrogen may be present without showing any spec- 
tral signatures and would serve to reduce the spectro- 
scopic temperature. In fact, hydrogen would go unno- 
ticed spectroscopically until at least log (H/He) = — 3.5 
was included. If we assume a log (H/Hc) = —4.5 (i.e., 
not spectrally visible), we obtain a spectroscopic solu- 
tion with Teff = 25,800 K, slightly more consistent with 
the photometric temperature. It is possible that trace 
amounts of hydrogen are present, however, without ad- 
ditional information (i.e., a trigonometric parallax), we 
have chosen to adopt the spectroscopic result assuming 
a pure helium composition. 

We note that the photometric temperature is far more 
uncertain. The optical photometry used to fit the SED 
is rather insensitive to effective temperature because the 
magnitudes fall largely on the Rayleigh- Jeans tail of the 
SED for such a hot object. Therefore, the spectroscop- 
ically determined effective temperature is significantly 
better constrained (spectrum and fit are plotted in Fig- 
ure 7). 

WD 2133-135: A DA WD that was imcovered dur- 
ing our sift of the SSS database. It is also included in 
Pauli ct al. (2006) (labeled HE 2133-1332 in that work), 
in which they analyzed spectra of 398 DA WDs to de- 
termine kinematic properties of the sample as part of 
the SN la Progenitor surveY (SPY). The authors state 
that their sample was generated from databases of known 
WDs; however, we performed a thorough search for a pre- 
vious spectroscopic confirmation of this object's WD sta- 
tus and were unsuccessful^. We have confirmed that this 
object was selected to be observed by the SPY project 
using data taken from the Hamburg/ESO (HE) survey 

^ Resources searched include the WD database of McCook 
& Sion (1999) (the current Web-based catalog can be found at 
http:/ /heasarc. nasa.gov/W3Browse/aII/mcksion. html) as well as a 
previous WD publication using data from the Hamburg/ESO (HE) 
survey (i.e., Christlieb et al. 2001). 



in which there was no previous confirmation spectrum 
but rather a high confidence in this object's luminos- 
ity class (R. Napiwotzki 2008, private communication). 
Thus, Pauli et al. (2006) is the first publication to con- 
firm this object's WD status (hence, it is a member of the 
known sample) yet this publication is the first to present 
a confirmation spectrum (see Figure 3). The authors 
estimate a spectroscopic distance of 23.3 pc, consistent 
with our distance estimate of 20.4 ± 3.5 pc. 

WD 2159-754: A known DA WD that has been ana- 
lyzed spectroscopically by Kawka et al. (2007) for which 
they determine the spectroscopic Teff = 9040 ± 80 K, 
moderately consistent with our photometric Teff = 9556 
± 242 K. In addition, they conclude that this object has 
a large surface gravity (log g = 8.95 ± 0.12) and hence 
a large mass (1.17 ± 0.04 Mq). Thus, their estimated 
distance (14 pc) is significantly closer than our estimated 
distance (30.5 ± 5.3 pc) based on an assumed log g = 
8.0. Trigonometric parallax measurements are underway 
to confirm its luminosity and hence its mass. 

5. DISCUSSION 

We continue to fill in the nearby WD sample with the 
discovery of 21 new WD systems brighter than V ^ 17. 
While the sample size is smaller than the number of new 
WDs presented in Paper XIX, the number of new sys- 
tems estimated to be within 25 pc is larger (eight re- 
ported in Paper XIX vs. twelve reported here for a total 
of 20 new WDs), due largely to the revised criteria set 
for candidates to be targeted for follow-up spectroscopic 
observations (see § 2). In addition, of the known samples 
evaluated in this effort, a total of 18 WDs are estimated 
to be within 25 pc (twelve reported in Paper XIX and 
six reported here). Combining the new and known sam- 
ples from both publications, we have found a total of 38 
WDs estimated to be within 25 pc - a volume that cur- 
rently contains 110 WD systems disregarding any quality 
constraints on the parallaxes.® Trigonometric parallax 
determinations are underway to confirm membership to 
the 25 pc sample but should the distance estimates prove 
accurate, the local sample of WDs will be increased by 
34%. Once these nearby WDs are confirmed, focused 
efforts can yield crucial WD masses and permit compan- 
ionship assessments. 

To ensure a reliable sample, we have begun to adopt 
the quality constraint that the trigonometric parallax er- 
ror cannot be greater than 10% of the parallax. Given 
the precision of ground-based trigonometric parallaxes 
(^2.0 mas or better), this constraint is entirely reason- 
able (even for a system at the 25 pc horizon, the 10% 
constraint amounts to an error of 4.0 mas). With the 
constraint in place, thirteen systems (five in the north 
and eight in the south) are eliminated from the known 25 
pc WD sample (now with a total of 97 systems). To bet- 
ter constrain their distances, we are currently measuring 

® This sample of WDs was compiled using parallaxes, both for 
the WDs as well as for additional components if the WD is part of a 
multiple system, from the Yale Parallax Catalog (van Altena ct al. 
1995), Hipparcos (Peiiyman et al. 1997; van Leeuwen 2007), and 
other recent efforts involving trigonometric parallaxes (e.g.. Smart 
et al. 2003; Gould & Chaname 2004; Ducourant et al. 2007). A 
comprehensive list, including weighted mean parallaxes when mul- 
tiple parallaxes arc avaible for a system, is included as an elec- 
tronic supplement to this publication and can also be found at 
http: / /www.DenseProject.com. 
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trigonometric parallaxes for seven of these systems in the 
southern hemisphere (one system, WD 1043—188 is too 
close to a companion whose brightness exceeds the bright 
limit of the 0.9m at CTIO). We would encourage north- 
ern hemisphere parallax programs to do the same for the 
five systems in the north (WD 0644+025, WD 0955+247, 
WD 1309+853, WD 1919+145, and WD 2117+539). 

If we separate new and known samples from Paper XIX 
and this work by proper motion, we see that the major- 
ity of added WD neighbors estimated to be within 25 
pc is found at lower proper motions. Table 2 breaks 
down this complete sample into proper motion bins of 
0.2" yr^^ below 1.0" yr""'^ (the first value in each col- 
umn corresponds to the systems presented in Paper XIX 
while the second value corresponds to those systems re- 
ported here). A quick summation shows the number of 
systems with /i > 0.6" yr^^ (12) is dwarfed by the num- 
ber of systems with /i < 0.6" yr~^ (26), including five 
that have fi < 0.2" yr~^ (of which, four have /x between 
0.18-0.2" yr~^). These results support the notion that a 
significant number of nearby WDs may still be found at 
very low proper motions. 

Other recent efforts have focused on the local WD pop- 
ulation. In particular, Holberg et al. (2008) have targeted 
the 20 pc sample determined by trigonometric parallaxes 
as well as photometric/spectroscopic parallaxes. The au- 
thors estimate the local WD density based on the 13 pc 
WD sample using the assumption that the 13 pc sample 
is largely complete. Paper XIX contained two objects 
(WD 0821-669 and WD 1202-232) previously unknown 
to be nearby, whose distance estimates as well as unpub- 
lished trigonometric parallaxes (Subasavage et al. 2008, 
in preparation) placed them within 13 pc. These have 
been included by Holberg et al. (2008) in their updated 
estimate of the local WD density. Given that none of 
the new WD discoveries in this publication have distance 
estimates within 13 pc, the notion that the 13 pc WD 
sample is largely complete is supported. The sample of 
Holberg et al. (2008) includes five new WD discoveries 
from Paper XIX whose distance estimates lie within 20 
pc (WD 0121-429, WD 0344+014, WD 0821-669, WD 
2008-600, WD 2138-332) as well as one from this pa- 
per (WD 0751-252, see § 4.2). The authors estimate 
that the 20 pc sample is ^80 % complete and that ^33 
± 13 WDs remain to be discovered between 13 and 20 
pc (again, assuming the 13 pc sample is complete). Five 
syst(;nis reported here (not including WD 0751—252) are 
estimated to lie within 20 pc. Thus, we are filling the 
incompleteness gap. 

A star's mass is one of its most important characteris- 
tics that, when coupled with composition and luminosity, 
defines several fundamental properties of that star such 
as its internal structure, future evolution, and total life- 
time. The same is true for WDs; however, only four WDs 
have measured astrometric mass determinations better 
than 5% (Sirius B, Procyon B, 40 Eri B, and Stein 2051 
B, Provencal et al. 1998). The identification of WDs in 
binaries, in particular, double degenerate systems, may 
increase the number of accurate astrometric masses of 
WDs if the system were resolvable using high-precision 
astrometric techniques (e.g., speckle, adaptive optics, or 
interferometry via Hubble Space Telescope's Fine Guid- 
ance Sensors). Of course, the nearer the system is to the 
Sun, the greater its projected separation and more likely 



TABLE 2 

Distance Estimate Statistics for New and Known White 
Dwarf Systems". 
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.. 1 + 
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The first number is from Paper XIX, and the second number is 
from this paper 



it is to be resolved, all else being equal. A sizable sam- 
ple of accurate WD masses will help constrain and revise 
WD models, while comprehensive searches for compan- 
ions in a volume-limited sample will allow an accurate 
multiplicity fraction to be determined. 

We continue to measure trigonometric parallaxes with 
good precision from the ground (~1.1 milliarcsecond er- 
rors on average, Subasavage et al. 2008, in preparation) 
for the nearby WD sample as part of the CTIOPI pro- 
gram. In addition, nearly all WDs within 15 pc in the 
southern hemisphere (including those we have presented 
in Paper XIX and here) have been targeted for long-term 
astrometric monitoring in search of perturbations from 
unseen companions. This effort is an extension of the 
Astrometric Search for Planets Encircling Nearby Stars 
(ASPENS, Koerner et al. 2003) that targets all red and 
white dwarfs within 10 pc in the southern hemisphere. 
Because of the spectral signatures of WDs (broad absorp- 
tion lines or no lines at all), astrometry is currently the 
only practical method for detecting sub-stellar/planetary 
companions to WDs. Unlike radial velocity variations 
used to detect planetary systems, astrometric signatures 
are linearly related to the distance to the system (the far- 
ther the system, the smaller the astrometric signature). 
Thus, a careful evaluation of the nearest WDs provides 
the highest probability for detecting astrometric signa- 
tures of companions to WDs. 
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7. APPENDIX 

In order to ensure correct cross-referencing of names for 
the new and known WD systems presented here, Table 
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table 3 astrometry and alternate designations for new and 
Known White Dwarfs. 



WD Name 


RA 


Dec 


PM 


PA 


Ref 


Alternate Names 






(J2000.0) 


(J2000.0) 


(arcsec yr ) 


(dcg) 






New Spectroscopically Confirmed White Dwarfs 


0011- 


-721 .... 


00 13 49.91 


-71 49 54.2 


0.326 


141.3 


S 


NLTT 681, LP 50-73 


0102- 


-579 .... 


01 04 12.14 


-57 42 48.6 


0.239 


091.1 


S 


SCR 0104-5742 


0123- 


-460 .... 


01 25 18.03 


-45 45 31.1 


0.759 


137.8 


s 


SCR 0125-4545 


0134- 


-177 .... 


01 37 15.16 


-17 27 22.6 


0.319 


189.3 


s 


NLTT 5424, LP 768-192 


0149- 


-723 .... 


01 50 38.49 


-72 07 16.7 


0.334 


223.9 


s 


SCR 0150-7207 


0311- 


-649 .... 


03 12 25.68 


-64 44 10.8 


0.190 


105.6 


s 


WT 106, LEHPM 1-3159 




— ODU .... 


riA QO 1^1=; 87 


QC C7 OQ Q 


U.oUl 




o 
o 




0431- 


-279 .... 


04 33 33.58 


-27 53 24.8 


0.403 


092.4 


s 


NLTT 13532, LP 890-39, LEHPM 2-405 


0620- 


-402 .... 


06 21 41.64 


-40 16 18.7 


0.379 


166.0 


s 


LEHPM 2-505 


0651- 


-398A... 


06 53 30.21 


-39 54 29.1 


0.227 


344.4 


s 


WT 202 


0651- 


-398B . . . 


06 53 35.34 


-39 55 33.3 


0.229 


344.7 


s 


WT 201 


0655- 


-390 .... 


06 57 05.90 


-39 09 35.7 


0.340 


242.6 


s 


NLTT 17220, L 454-9, LTT 2692 


0708- 


-670 .... 


07 08 52.28 


-67 06 31.4 


0.246 


246.3 


s 


SCR 0708-6706 


0709- 


-252 .... 


07 11 14.39 


-25 18 15.0 


0.223 


334.4 


s 


SCR 0711-2518 


0751- 


-252 .... 


07 53 56.61 


-25 24 01.4 


0.426 


300.2 


s 


SCR 0753-2524 


0816- 


-310 .... 


08 18 40.26 


-31 10 20.3 


0.842 


162.6 


s 


SCR 0818-3110 


0856- 


-007 .... 


08 59 12.91 


-00 58 42.9 


0.202 


125.8 


s 


NLTT 20690, LP 606-32 


1116- 


-470 .... 


11 18 27.20 


-47 21 57.0 


0.322 


275.1 


s 


SCR 1118-4721 


1817- 


-598 .... 


18 21 59.54 


-59 51 48.5 


0.365 


194.9 


s 


SCR 1821-5951 


1916- 


-362 .... 


19 20 02.83 


-36 11 02.7 


0.208 


132.0 


s 


SCR 1920-3611 


2009- 


-471 .... 


20 12 48.75 


-46 59 02.5 


0.244 


136.3 


s 


WT 689 


2118- 


-388 .... 


21 22 05.59 


-38 38 34.7 


0.186 


113.5 


s 


SCR 2122-3838 


Known White Dwarfs (without Trigonometric Parallaxes) 


0233- 


-242 .... 


02 35 21.80 


-24 00 47.3 


0.620 


189.5 


s 


LHS 1421, NLTT 8435, LP 830-14 


0707- 


-320 .... 


07 09 25.07 


-32 05 07.3 


0.551 


338.2 


L 


LHS 1898, NLTT 17486, LP 896-18 


1223- 


-659 .... 


12 26 42.02 


-66 12 18.5 


0.190 


182.0 


L 


NLTT 30737, LP 104-2, GJ 2092 


1241- 


-798 .... 


12 44 52.70 


-80 09 27.8 


0.578 


306.3 


L 


LHS 2621, NLTT 31694, LP 38-80 


2007- 


-219 .... 


20 10 17.51 


-21 46 45.6 


0.311 


158.0 


L 


NLTT 48815, LP 870-43 


2133- 


-135 .... 


21 36 16.38 


-13 18 34.5 


0.297 


120.2 


s 


NLTT 51636, Ross 203, HE 2133-1332 


2159- 


-754 .... 


22 04 20.84 


-75 13 26.1 


0.523 


277.6 


s 


LHS 3752, NLTT 52728, LP 48-15 


2216- 


-657 .... 


22 19 48.31 


-65 29 17.6 


0.660 


160.1 


L 


LHS 3794, NLTT 53489, LP 119-34 


2306- 


-220 .... 


23 08 40.78 


-21 44 59.6 


0.350 


109.0 


s 


NLTT 55932, LP 877-69 


2336- 


-079 .... 


23 38 50.74 


-07 41 19.9 


0.034 


126.6 


Su 


GD 1212 


2351- 


— o.'"S5 .... 


2:-! 51 01.11 


-:-!:-! 1() :-10.:-! 


0.500 


21(),5 


L 
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